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Abstract: 

We present a two channel eikonal calculation in which the rescattering through a diffrac- 
tive channel is included in addition to the elastic channel. Considering the spread of the 
experimental data, we find that we can obtain a very good description of a tat , a e i and B d in 
the ISR - Tevatron energy range. In this range of energy the diflractive channel, that was 
included in our calculation, leads to a ratio of o snivel which varies between 1 and 0.5 for 
20 GeV < \/s < 14 TeV in agreement with the experimental data. The calculated survival 
probability of dijet production with a large rapidity gap is consistent with the data. 
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1 Introduction 



A simple Regge-pole parameterization suggested by Donnachie and Landshoff (DL)[l|], has 
achieved remarkable success in describing hadronic and photon induced total and elastic 
cross sections. The elastic amplitude is, obviously, bounded by s-channel unitarity requiring 
screening (shadowing) corrections (SC) at sufficiently high energies. However, up to the 
Tevatron energies no decisive experimental signature of this effect has, thus far, been seen 
in the elastic channel. 

An overall view of the features of high energy collisions @ || [0 f| || reveals at least 
two channels which are not compatible with the DL parameterization. 



1. Since the total cross section behaves at high energies as s Ap , where «p(0) = 1 + Ap, 
we expect both a e i and usd, the integrated single diffraction (SD) cross section, to 
behave approximately like s 2Ap (since the Pomeron trajectory slope is small we can 
neglect the effects of the t integration). This is, indeed, the experimentally observed 
behaviour of a e i, but ctsd has a much milder energy dependence clearly seen in the 
ISR- Tevatron energy range 0. 

2. < \S\ 2 >, which is definedf| as the survival probability of two jets with a large rapidity 
gap (LRG), is measured [f|]|f| to be a small fraction of unity decreasing with energy. In 
a DL type model, containing very weak rescattering corrections, we expect < |5| 2 > 
to be approximately one with no, or very little, dependence on energy. 



These problems have been discussed in our previous piiblicatioiis0 (jlTj ]TT[ ]T2[] JT^] where 
we have utilized a simplified one channel eikonal approximation in which the b-space scat- 
tering amplitude is given by 



a(s, b) 



n(s,b) 



the opacity is assumed to have a Gaussian form 



n(s,b) 



u(s)e . 



(2) 



Even though this may be considered as an over simplified toy model, it provides a semi 
realistic reproduction of the main experimental features. Our main conclusion [|TTJ fl2| is that 
while SC saturate the elastic channel in the TeV range, diffractive channels are saturated 
at relatively low (ISR) energies. This feature results in a mild energy dependence of asn at 
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fixed M 2 , where the power behaviour of the cross section with s approaches a Ins dependence 
in the high energy limit. The above eikonal approximation also leads |]1U||13|| to estimates of 
< | S | 2 > which are compatible with the Tevatron[|J and HERA|| data. 

The one channel eikonal approach contains two severe deficiences:: 



1. The main input assumption in such a model is that the rescattering is elastic, i.e. 
a diff « °ei- This is clearly not corroborated by the data in the ISR-Tevatron energy 
range. 

2. The introduction of SC introduces a fundamental problem in the definition of <Jsd- in 
a missing mass experiment, the integrated SD cross section is defined as 

rt max (M 2 ) rM^ ax d 2 <7 S D 

11 SD Jt mi n(M*) JMl in dM 2 dt ' 

where the upper M 2 integration limit is taken as a fixed fraction of s, usually 0.05s. 
In a non screened triple Regge approximation JTJj <jsd behaves like s 2ap ^~ 2 at fixed 
M 2 . With a super critical Pomeron (a p > 1), the M 2 integration converges, leading 
to an approximate s 2Ap behaviour of the integrated cross section. However, once SC 
are introduced [|ITJ, the M 2 dependence of cjsd changes in the asymptotic limit to Ins. 
Accordingly, the M 2 integration is divergent in s due to the upper M 2 integration 
limit. The predictions are, thus, flawed if examined over a wide energy range. Indeed, 
while we managed to fit gsd in the UA(4)-Tevatron range, our calculated a$D values 



lTJ in the ISR energy range fall below the data 



2 The Model 

Our goal is to construct an eikonal model in which the rescattering can be either elastic or 
diffractive. In the simplest approximation we consider the diffractively produced hadrons as 
a single hadronic state. We have, thus, two orthogonal wave functions 

< V h \V D > = (4) 

where ^ is the wave function of the incoming hadron, and *S?d is the wave function of the 
outgoing diffractively produced particles, initiated by this hadron. 

Consider two wave functions and ^ which are diagonal with respect to T, the inter- 
action operator. The amplitude of the high energy interaction is equal to 

A i>k = < tfitffclTl^/tftf > = A i:k 5 iti ,5 kr . (5) 
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In a two channel model i, k = 1, 2. The extension to a higher number of channels is straight 
forward. The amplitudes A i>k satisfy the diagonal unitarity condition 

2ImAf >k (s,b) = \Af >k (s,b)\ 2 + G7? k (s,b) (6) 
for which we write the solution 

A( k {8,b) = i (l - e-^J ; (7) 
G&(s,&t) = 1 - e -n*.*(».*), (8) 
where Jlj^ is the opacity of the (i, fc)-th channel with a wave function x ty k . 
In this representation ^ and typ can be written as 

^ h = a*i + /9^ 2 ; (9) 
\£ D = + (10) 

Since l^/J 2 = 1, we have 

a 2 + f3 2 = 1 (11) 
The wave function of the final state is 

% = \T\V h x V h > = (12) 

a 2 ^!^! X + a/3A li2 {^i x # 2 + tf 2 x #i} + /3 2 A 2;2 *2 x ^ 2 
We now define 

a eZ (s,6) = < * h x > = a 4 + 2a 2 / 3 2 A lj2 + P A A 2>2 ; (13) 

osd(s,&) = < x > = aP{ -a 2 A hl + (a 2 - /3 2 ) A 1>2 + /3 2 A 2)2 } (14) 

odd(s,&) = < x * D |*j > =a 2 (3 2 { A 1A - 2A lj2 + A 2j2 } . (15) 

Assuming, that the double diffract ive production is small, we have 

a el (s,b) = A 1A - 2f3 2 (A M - A 1<2 ) ; (16) 
a SD {s,b) = -af3(A 1A - A lj2 ) . (17) 

a e i and asD can be written in terms of the opacities Qi = Q^i and Q 2 = Jli ;2 

a ei (s,6) = i ( 1 - e 2)J-2^ 2 [e 2 - e 2 ) ; ( 18 ) 
a S D{s,b) = -a(3 \e 2 - e 2 I. (19) 



Defining AQ = fi 2 — fh, we get 

a e i(s,b) = ill-e 2 )j-2/5 2 e 2 I 1 - e 2 ; (20) 
asD{s,b) = -a fie 2 ( 1 ~ e 2 J • ( 21 ) 

In the limit where /3 << 1 and Af2 << 1, we reproduce the single channel eikonal model 

(ni (3,6) \ 
1 - e — )J ; (22) 

= (3^e-^. (23) 

For a convenient semi realistic parameterization, we follow our earlier publications and 
adopt a Gaussian approximation for the opacities 

nf( a ,6) = ^ (i-)^ e -(W = ^(-) A ^"^W; (24) 

An p (s,b) = -Jf „ 2 (-) A ^e'^ = <(-) A -e~^\ (25) 
For both radii we assume the form 

(Rff = «? + 4a£, ln( S / So ) , (26) 
where % = 1,D, and take s = 1 GeV 2 . 

Since our investigation covers the ISR energy range , we need to include, in addition to 
the Pomeron parameters defined above, a Regge contribution which is defined in a similar 

way 



nf(a,6) = ^ = ^(-)-"e ^«> a ; (27) 

n-B o b2 t£ 



AQR ^ = ^T„ 2 (-r^ (flgW)2 = ««S«' a ; (28) 

1T{K D {S)) S S 

{Rff = «) 2 + ln(s/ So ) ; (29) 

where 77 = 1 — or(0). 
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In the following we check the ability of this simple model to describe the high energy 
data on 



Otot(s) = 



2nJ™db 2 a el (s,b); 
tt / °° db 2 \a el (s,b)\ 2 : 




(30) 
(31) 

(32) 

(33) 



B el (s) 
&sd(s) = 



/ Q °°db 2 b 2 a el (s,b) _ 



2 db 2 a el (s,b) ' 

2-k ^db 2 \a SD {sM 



where we use Qi(s, b) = Qf(s,b) + f2f(s,6). 

We assume single diffraction to be the only non negligible diffractive channel. By av- 
eraging the pp and pp data, we eliminate the need to discuss odd parity Regge exchanges, 
which may be marginally important at the lower ISR energies. The elimination of the odd 
contribution is exact for a tot and approximate for the elastic channel. 

Once the parameters of our model are specified we can calculate the survival probability 
[BJ of dijets with a LRG. We consider this is an important consistency check and discuss it 
in Section 4. 



3 Data base and fitting procedures 

Our study aims at a better understanding of the role s-channel unitarity screening and its 
relevance to soft Pomeron physics. Accordingly, we have limited our investigation to high 
energy pp and pp data at the ISR and above, i.e. s > 300 GeV 2 . This choice is made so as 
to minimize the dependence of our analysis on the secondary Regge exchanges. 

Our data base contains 61 entries: 

1. 18 values for o to t [@]- In the ISR range we have 10 points averaged between pp and pp. 
Above ISR there are 8 measured values of a to t{pp). These energies are sufficiently high 
to neglect the odd parity contributions. 

2. 9 a e i data points[|[, 5 of which are averaged ISR values and 4 higher energy cross 



sections. 
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11 values for 5 e /@. For 2 ISR energies, where we have the data, we took the averaged 
value. 
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21 <jsd data points^, all of which were obtained from missing mass experiments 
initiated by pp. The data corresponds to the sum of SD produced at either the proton 
or anti proton vertex. We have omitted from the analysis the SD data reported by 
Ref. ||15||. These data points are systematically much lower than the other ISR reported 



SD cross sections [y]. 

Our model, as specified in Section 2, contains Pomeron and effective Regge trajectories 

a P (t) = 1 + A P + 0.2 1 ; (34) 

a R (t) = 0.6 + t ; (35) 

where Ap is a fitted parameter. The assumed values of a p , or(0) and a R are those commonly 
used in Regge phenomenology. The other fitted parameters, defined in Section 2, are (3, Ap, 

a oD- I n accordance with the triple Regge formalism we take 

(Rl D ) 2 = ^oi) 2 where i = P, R, (36) 



We neglect the triple vertex radius [p[ 



The predictions of our best fit, compared with the data are presented in Figs 1-4. The 
values of the fitted parameters are given in Table I. Our best fit has an overall x 2 = 129.8 

2 

corresponding to -jj- = 2.6. This is a non satisfactory high value. However, when inspecting 
the x 2 contribution of the individual points we find that 7 points produce x 2 = 65.2. Ne- 
glecting these 7 points results in -fj = 1.5. Re-fitting the data without these 7 points, we 
get essentially the same fit. Conidering the spread of the experimental points reported with 
small quoted errors, we consider our fit to be a good one. 



4 Survival probabilities 

A LRG process is identified by the absence of produced hadrons in a sufficiently large rapidity 
gap region. This is regarded as a reliable signature for a crossed channel exchange of a 
Pomeron or alternatively a colourless gluonic state, be it non perturbative or perturbative. 
Recent experimental investigations of jets produced in the Tevatron(|]f^] and HERA|J0 
have recorded this phenomena. We shall discuss here the simplest process where we have 
two jets produced with a LRG separation. 

The experimentally measured quantity of interest is f gap , which is the ratio of the cross 
section of dijet production with a LRG, and the inclusive two jet cross section (see Fig. 5) 

f 9 a P = ° lM2 36tS ] = <| S | 2 > -F s , (37) 
Jgap a md (2 jets) \ \ 
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Figure 3: Slope of the elastic cross section versus energy. 
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Parameters 


Best fit values 


A P 


0.126 


P 


0.464 




12.99 


(R^(GeV- 2 ) 


16.34 


< D ( GeV^ ) 


145.6 




4.78 


(Roi) 2 (GeV~ 2 ) 


12.44 


a* D ( GeV- 2 ) 


999.0 



Table 1: The values of the param- 
eters of the Pomeron and Reggeon 
input Eikonal amplitudes. 
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Figure 5: Pictorial definition of f gap , where P and G represent, respectively, the exchange of 
a colour singlet and a colour octet. 

F s is the dynamical ratio of singlet to octet colour exchange leading to two jet production 
and its calculation is beyond the scope of this paper. F s has to be modified by the survival 
probability < \S\ 2 >, which is the probability that the produced LRG event is not filled by 
hadronic debris resulting from rescattering of partons or hadrons. 

The calculation of < \S\ 2 > is model dependent. We shall follow the procedure suggested 
by BjorkenQ, in which the survival probability is the normalized b-space convolution of 
the hard partonic scattering amplitude and P(s, b) - the probability that the two initial 
projectiles do not interact inelastically. 



Taking a Gaussian approximation for both Q(s, b) and a#(s, b) simplifies the calculation and 
we obtain 




(38) 



In the following we neglect the dependence of a# on the rapidity gap of interest. 
In a single channel eikonal model 



P(s,b) = e- n{s ' b) . 



(39) 




(40) 
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where j(a, x) = Jq z a 1 e z dz, v is defined in Eq. (§) and a = jgf ^ . 

The results of our single channel eikonal calculation [[HJ are in fair agreement with the 
DO datajl. However, the calculated energy dependence is not sufficiently strong. This 



deficiency is removed once we estimate [|l^] u, R 2 and R 2 H directly from the data. The two 
channel eikonal model presented in this paper may have some problems in reproducing the 
experimental energy dependence, since the data for <Te ' +CT ' M// is almost energy independent 
in the ISR-Tevatron energy range, whereas ^ grows monotonically with s. For this reason 
we consider the calculation of < \S\ 2 > in a multi channel eikonal model, to be an important 
check of the validity of this approach to calculating the SC. 

A hard LRG b-space amplitude, uncorrected by the survival probability, is given in our 
formalism by 

a H = a 2 a? + f3 2 af (41) 

where af^ and at? correspond to processes initiated by our and ^2 wave functions. As 
previously, we impose the restriction that no inelastic rescattering modifies our hard process, 
by multiplying each component of an by the corresponding e~ Qn . 

The survival probability in our two channel model is given by: 

<| S ( S )|>>= (42) 



where 



7T 

and 



D(s) = vr jH db 2 \ Y, a? P t (s,b)j = (43) 
| o °° db 2 (e-W){l - (3 2 [l - e- An{s ^]}af + (3 2 e~ Q ^ e~ An ^ af D ) ; 

f'OO { ^ \ c 00 

n(s) = 7T y o db 2 r£ of J = 7r j / o db2 { a i + a s D ) ■ ( 44 ) 

<| S{ s ) \ 2 > can be calculated numerically once our free parameters have been fixed 



provided we know R 2 H . We have checked that taking R 2 H = 8 GeV -2 [|16|, we obtain <| 
S( s) \ 2 > values in the Tevatron range which are in a fair agreement with the DO data [|J. 
The calculated energy dependence of <| S( s) | 2 > is given by 

1.3 - 1.4 

to be compared with DO value of 2.2 ± 0.8. In as much as this is just a consistency check 
and not an actual fit, we consider these results to be reasonable. A improved and detailed 



< 


1 S( y/s = 630 GeV) | 


2 > 


< 


S(y/s = 1800 GeV) 


\ 2 > 



analysis of LRG events in a multi channel eikonal model will be published soon [17| • 
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5 Discussion 



The following are the main features and predictions of our suggested model. 



1. The consistent treatment of the diffractive dissociation channel in our model leads to 
a good description of the experimental data on the single diffraction in ISR - Tevatron 
energy range. We reproduce the energy dependence of the ratio o sol^ei as well as its 
value ( see Fig. 8a ) in agreement with the experimental data. 

2. Our input Pomeron has Ap = 0.126 which determines our high energy predictions 

( see Figs. 1 - 4 ). For LHC, yfl = 14TeV, we predict a tot = 104.3 mb. Our 
prediction is slightly higher than the DL prediction of 101.5 mb and is compatible 
with a recent prediction of Block et al.||18|| of 108 ± 3.4 mb. The model suggested in 
Ref. [0 is a single channel eikonal calculation with a Regge input. DL approach 
is a non screened Regge model corrected by a weak Pomeron - Pomeron cut. The 
small differences between the a to t predictions at LHC energies of these rather different 
models, re-enforces our observation that SC are weak for a to t- 
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20.5 GeV~ 2 . 
This is com- 
patible with the small difference between the two a to t predictions, combined with our 
predicted value of B e \ being higher than the prediction of Block et al. 



Our LHC predictions for the elastic channel are a e i = 24.5 mb and B e \ 
Our a e i prediction is about 15% lower than the prediction of Ref. 



4. Some insight into the difference between a non screened DL type 6-space elastic am- 
plitude, and our screened a e i(s,b) is provided Fig. 6. As seen in Fig. 6a the DL non 
screened amplitude violates s-channel unitarity ( to be distinguished from the Frois- 
sart bound ) at small impact parameters starting from yfs ~ 3 TeV. The introduction 
of a weak P-P cut [[!]] reduces a^ L (s, b) at small b by approximately 10%, leaving the 
DL prediction for a tot practically unchanged. In comparison, our screened a el (s,b) is 
considerably lower and wider ( see Fig. 6b ). This behaviour is a consequence of both 
screening and the existence of a competing diffractive channel. 



5. As we have noted, the investigation of <jsd in the ISR- Tevatron energy range provides 
strong support for the importance of SC at existing collider energies. It is interesting to 
compare our approach with the Pomeron flux renormalization suggested by Goulianos 
[ 19| . In the eikonal model, the damping due to unitarity in the diffractive channel is 
given by e~ n *^ 3,b \ where i denotes the rescattering channel considered. As such, the 
damping is 6-dependent, with maximum suppression at small b resulting from the fact 
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that the input Qi(s,b) are central. As we saw, the damping is smooth and becomes 



significant at ISR energies. In Ref.|T3, a unitarity flux correction is applied only when 
the Pomeron flux exceeds unity. This amelioration is b- independent. Numerically, our 
(Tsd has a stronger energy dependence than suggested in Ref. [JHJ. In the ISR-Tevatron 
range the difference between the predictions is small. However, at LHC, \/s = 14TeV, 
we predict <jsd = 12 mb to be compared with Goulianos' estimate of 10.5 mb. 

A simple Gaussian input for the opacity leads to dips in at lower values of t than in 
the data. To obtain these dips at their observed experimental values, the Gaussian form 
of Q(s, b) has to be replaced by a dipole formf2(|. We consider this to be a relatively 



minor technical modification of the calculation. We note that the dip structure can be 
reproduced in an eikonal calculation |2(J, as well as by the introduction of a weak P-P 
cutn, leaving the DL results on o tot and o e \ virtually unchanged. 

The distinctive features of our diffractive 6-space amplitude are visible in Fig. 7. Our 



g sd( s 5 b) is peripheral. This is to be compared with Goulianos et al. who produce 
a SD amplitude proportional to the unscreened DL amplitude shown in Fig. 6a. Our 
model is, thus, compatible with the Pumplin's bound []2T] 



\a D (s,b)\ 2 = \a el (s,b)\ 2 + \a SD (s,b)\ 2 < 1, (45) 
which translates, after a 6-integration, to 

a d (s) + cr SD (s) < ^cr tot (s) . (46) 

In Fig. 8 we show the energy dependence of R e i = - £eL , Rsd = £j2s and Rd = R e i + 
Rsd ( Fig.8a ), and the energy dependence |azj(s, 6) | 2 = \a e i(s, b)\ 2 + |asr>(s, &)| 2 
at b =0 ( Fig. 8b ). All ratios are well below the Pumplin bound. We have checked 
that the bound is saturated at the non realistic ultra high energies of ^Js ~ 10 9 GeV. 
The Goulianos model differs in as much as it violates the Pumplin bound at small b. 
This is a consequence of both a e i(s,b) and aso(s,b) peaking at b = 0. We note that 
the Pumplin's bound was proven in a multichannel eikonal model. Its validity in other 



screening models is not clear p2[ . A peripheral 6-space structure for asD implies dips 
or breaks in at small t, where the details depend on form of the input taken in the 
eikonal calculation. Some indications for such a structure have been observed in ISR 



SD experiment with a very small diffractive mass p3[ . More accurate experimental 
data are required to investigate this aspect further. 

8. A more detailed treatment of <| S \ 2 > in a multichannel eikonal model is in preparation 



171 . In the present paper we have calculated this observable as a check of our model. 
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We obtain a reasonable energy dependence of <| S | 2 > even though Rd is almost energy 
independent in the ISR- Tevatron energy range ( see Fig.8a ). The resulting energy 
dependence of <| S | 2 > in our model is not surprising as each channel is suppresed by 
e -n i (s,6)_ ^ g ] on g as gQ are small, the energy dependence of <| 5* | 2 > is determined 
by Rd which is almost flat in the energy range of interest. When the SC become 
significant, the relevant ratios are R e \ and Rsd, both of which show sufficient energy 
dependence to induce a final <| S | 2 > which is compatible with the DO measurements 



9. In a recent preprint Cudell et al. ||24|| come to a conclusion similar to ours for a tot- 
However, Cudell et al. use this result to question whether unitarity corrections are 
necessary at all. We wish to emphasize again that the value of SC differs for different 
channels. They are mild for the elastic channels but appreciable for <jsd an d <| S | 2 >. 
Since we believe that the same colour singlet ( Pomeron ) mediates all these reactions, 
we advocate not drawing general conclusions which are based on the finding in one 
channel only. 
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Fig. 6-a Fig. 6-b 



Figure 6: The b - dependence of the elastic amplitude for the unscreened DL Pomeron ( 
Fig. 6a ) and in our two channel model ( Fig. 6b ). 
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Figure 7: The b - dependence of the single diffraction amplitude in two channel model. 
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Figure 8: Energy behaviour of the ratios R D = (a e i + CsD)/< J tot, Rei = °~ei/°~tot and 
Rsd = 0-sd/o-m (Fig. 8a ) and \a D (s,b = 0)| 2 = \a el (s,b = 0)| 2 + \a SD (s,b = 0)| 2 ( Fig. 8b). 
The dotted line in both figures shows the Pumplin bound. 
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